TOTAL POSITIVITY OF A CAUCHY KERNEL 



THOMAS SIMON 

Abstract. We study the total positivity of the kernel l/(a;^ + 2cos(7rQ)xj/ + t/^). The case of infinite 
order is characterized by an application of Schoenberg's theorem. We then give necessary conditions 
for the cases of any given finite order with the help of Chebyshev polynomials of the second kind. 
Sufficient conditions for the finite order cases are also obtained, thanks to Propp's formula for the 
Izergin-Korepin determinant. As a by-product, we give a partial answer to a question of Karlin on 
positive stable semi-groups. 



1. Introduction 

Let I he some real interval and K some real kernel defined on / x /. The kernel K is called 
totally positive of order n (TP„) if 

det [K{xi,yj)]^^. .^^ > 

for every m G {1, . . . ,n}, xi < ... < Xm and yi < . . . < i/m- If these inequalities hold for all 
n one says that K is TPqo- The kernel K is called sign-regular of order n (SR^) if there exists 
{£m}i<m<n G {"l, 1} such that 

for every m G {1, . . . , n}, xi < . . . < Xm and yi < . . . < i/m- If these inequalities hold for all n one 
says that K is SRoo- The above four properties are called strict, with corresponding notations STP 
and SSR, when all involved inequalities are strict. We refer to [7] for the classic account on this 
field and its various connections with analysis, especially Descartes' rule of signs. We also mention 
the recent monograph [TT] for a more linear algebraic point of view and updated references. 

A function / : M — t- M"*" is called Polya frequency of order n < oo (PF„) if the kernel K{x, y) = 
f{x — y) is TP„ on M X M. When this kernel is STP„, we will use the notation / G SPF„. Probability 
densities belonging to the class PFqo have been characterized by Schoenberg - see e.g. Theorem 
7.3.2 (i) p. 345 in [7] - through the meromorphic extension of their Laplace transform, whose 
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reciprocal is an entire function having the following Hadamard factorization: 
(1) = e-^^^^'^n(l + «n.)e-- 

with X the associated random variable, 7 > 0, (5 G M, and < 00. The classical example is 

the Gaussian density. PF2 functions are easily characterized by the log-concavity on their support 
- see Theorems 4.1.8 and 4.1.9 in [7J. But except for the cases n = 2 or n = 00 there is no handy 
criterion for testing the PF„ character of a given function resp. the TP„ character of a given kernel, 
and such questions might be difficult. In this paper we consider the kernel 

Ka{x,y) = ^ 

+ z cos(7rajxy + 

over (0, +00) X (0, +00), with a G [0, 1). Because of its proximity with the standard Cauchy density 
we may call K^, a Cauchy kernel, eventhough the denomination Poisson kernel would also be 
justified. Occurrences of the kernel in the literature are many and varied. We show the 
following 

Theorem, (a) One has 

gSTPoo ^ gSRoo ^ aG{l/2,l/3,...,l/n,...,0}. 

(b) For every n > 2, one has 

Ka e TP„ Ka e SRn ^ o G { 1/2, 1/3, . . . , 1/n} or a< 1/n. 

(c) For every n > 2, one has q < 1/n A l/(n^ — n — 6)+ =^ Ka G STP„. 

The well-known fact that Kq and are STPqo is a direct consequence of explicit classical 

formulae for the involved determinants, due respectively to Cauchy and Borchardt - see e.g. (2.7) 
and (3.9) in [8j - and which will be recalled thereafter. The characterization obtained in Part (a) 
follows without difficulty from Schoenberg's theorem and Gauss' multiplication formula, and will 
be given in Section 2. The more involved proofs of Part (b) and Part (c) rely both on an analysis 
of the derivative determinant 

r d^~^^~^Ka ' 
" [dxi-^dy^-\ 

where, here and throughout, we set det„ for the determinant of a matrix whose rows and columns 
are indexed by (i, j) G {1, . . . , n}^. To obtain Part (b) we establish a closed expression for AJJ(1, 0+) 
in terms of Chebyshev polynomials of the second kind, an expression which is negative whenever 
a {1/2, 1/3, . . . , 1/n} or a > 1/n. The computations are performed in Section 3, in four manners 
involving respectively Wronskians, Schur functions, rectangular matrices and alternating sign ma- 
trices. The observation that these four very different approaches all lead to the same formula was 
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interesting to the author. On the other hand only the last approach, which is based on Propp's 
formula for the Izergin-Korepin determinant, seems successful to get the more difficult Part (c). 
This latter result is partial because our upper condition on a is probably not optimal. We raise the 

Conjecture. For every n> 2, one has 

(2) a <l/n ^ Ka e STP„. 

This would show that the inclusion in Part (b) is actually an equivalence, which is true for 
n = 2, 3 by Part (c) of the theorem because then n > (n? — n — 6)+. In Section 5 we show ([2]) 
for n = 4, 5, and we also give some heuristics reasons supporting the validity of this inclusion for 
all n. The latter seems however to require wild computations. In Section 5, we state three other 
conjectures of combinatorial nature whose fulfilment would entail ([2]). The most natural one is a 
criterion for the positivity of the generating function 

ti(A) = k 

evaluated at a certain complex number z, where An is the set of alternating sign matrices of size n, 
A'^ is the anticlockwise quarter-turn rotation of A, fJ.{A) is the number of negative entries and iy{A) 
the inversion number (see the precise notations below). Whereas we can show this criterion for 
k = or k = /imax (see below Propositions 3 and 4), unfortunately we cannot do so for all k since 
there is no sufficiently explicit general formula for fn,k- Let us stress that the single evaluation of 
/n,fc(l) = '^{A G An, fJ-iA) = A;} is a difficult open problem, solved only for certain values of A; - see 
|10j and the references therein. 

The present paper was initially motivated by a question of S. Karlin - see Section 6 below 
for its precise statement - on the total positivity in space-time of the positive stable semi-group 
{t, x) I—)- Pait, x) on (0, -|-oo) X (0, -|-oo), which we recall to be defined by the Laplace transform 

/•oo 

(3) / pa{t,x)e-^''dx = e-^^", A > 0. 
As a simple consequence of Part (b), in Section 6 we show the 

Corollary. For every n > 2, one has 

Pa G SRn =^ a e {1/2, 1/3, ... , 1/n} or a < 1/n. 

The inclusion a < 1/n =^ G TP„ was proved in [12j for n = 2 and we believe that it is 
true in general. This would give a complete answer to Karlin's question and also entail the above 
conjecture - see Section 6 for an explanation. This will be the matter of further research. 
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2. Proof of part (a) 

This easy part of the theorem is a consequence of the following Proposition 1 and Corollary 1. 



Proposition 1. One has 

Ka gSTPoo ^ GTPoo ^ QG{l/2,l/3,...,l/n,...,0}. 

Proof. We begin with the second equivalence. Consider the generalized logistic distribution with 
density 

sin(7rQ) 

° 27ra(cosh(x) + cos(7ra)) 

over M. Simple transformations - see Theorem 1.2.1. in \J\ - entail that the TPqo character of 

amounts to the fact that ga G PFoo- For every s G (—1,1), compute 

sinfvra) 1'°° sinfvras) 

/ — 7^ -, ; ^du — 

lo 



na Jq V? + 2ucos(7ra) + 1 asin(7rs) 
where the right-hand side follows from the residue theorem and is meant as a limit for a = 0. If 
a ^ {1/2, 1/3, ... , 1/n, . . . , 0}, the function 

a sinfvrs) 
s ^ ■ ' ' 
sin(7rasj 

has a pole at 1/a so that PFqo by the aforementioned Theorem 7.3.2 (i) in [7]. If a = 1/n for 
some n > 2, writing 

sin(7rs) ^ r(l-g)r(l + g) 
nsin(7rs/n) T[l — s/n)T{l + s/n) 

and applying Gauss' multiplication formula and Weierstrass formula for the Gamma function - see 

e.g. 1.2(11) p. 4 and 1.1(3) p.l in [5j - shows that this function is of the type ([T]), in other words 

that ga € PFqo. If a = 0, the same conclusion holds true thanks to the Eulerian formula 

/ o2 ■ 

n 1 



sm vrs 



n>l 



This finishes the proof of the second equivalence. To show the first one, it remains to prove that 
Ka G STPoo whenever a £ {1/2, 1/3, . . . , 1/n, . . . , 0}. Cauchy's double alternant formula (see e.g. 
(2.7) in [8] or Example 4.3 in [IT] ') 



det^ 



1 



ni<i<j<n(2/j yi){^'j ^i) 
Y[l<i,j<n(^i + 



entails immediately that K-^,^ G STP^o- Analogously, Borchardt's formula (see e.g. (3.9) in [5]) 



(5) 



det^ 



.{xi + yjY' 



det^ 



1 

Xi + yj 



X perm„ 



1 

Xi + yj 
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yields Kq G STP^o. An alternative way to prove these two latter facts consists in writing 

^ = / e"^ "e"^^ du and = / e-^^e""^ udu. 

X + r Jo \x + V) Jo 

Indeed, the composition formula - see Lemma 3.1.1 in [7] - entails that Kq and STPoo 

because the kernel e~"^ is SRRoo on (0, +oo) x (0, +oo) - see [7] p. 18 for the latter fact and p. 12 

for an explanation of the notation SRRqo. 

The remaining cases a {0,1/2} are slightly more involved. First, it follows from (jj]) and a 

fractional moment identification that /i/„ is the density of the independent sum 

(6) iog(ri/„) + • • • + iog(r(„_i)/„) - iog(ri/„) - • • • - iog(r(„_i)/„) 

where Tt is for every t > the random variable with density 

1 g X 

Second, it is straightforward that the density of log(rf) is SPFqo. This property conveys to 
thanks to the above factorization and the composition formula. 

□ 

Remark 1. (a) The above argument based on the composition formula yields the STPqo character 
of all kernels {ax + by + c)~'^ for a,b,d > and c > (this is the main result of [4] - see Theorem 
3.1 therein), in writing 

[ax + by + c)'^ Jq 

(b) It is easy to see that x i— )• Ka{\/x, y^) is a completely monotone function for every y > 0, and 
that there exists a certain positive finite kernel La{x,y) on (0,+oo) x (0,+oo) such that 

POD 

Ka{x,y) = / e''^^'^ La{u,y)du. 



Jo 

When a is the reciprocal of an integer, by ([6]) it is possible to express Lq, as a convolution of 
weighted Laplace transformation kernels and show that it is RRoo- The kernel is less explicit in 
the other cases but we feel that its sign-regularity index matches the total positivity index of Kq,, 
in other terms that 

Ka € TP„ La £ RRn 

for all n. 



Proposition 2. The function is positive- definite for every a G [0, 1). 
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Proof. The beginning of the proof of Proposition 1 entails by analytic continuation that 

i,^ , , , sinh(7ras) 
asinh(7rsj 

for every s G M. Since a S [0, 1), we can apply the Fourier inversion formula to obtain 

sinhfvras) , ^ 



27rasinh(7rs 

The conclusion follows from Bochner's theorem. 



□ 



Remark 2. The above proposition is very well-known, but we gave a proof for the reader's comfort. 
It is also true that (7^ is positive-definite for every a G [0, 1) and r > - see e.g. Exercise 5.6.22 (ii) 
in [2]. 

Corollary 1. For every a € [0, 1) and n > 2, one has 

Ka e TP„ ^ Ka e SRn. 

Proof. Proposition 2 entails that for every a S [0, 1) the function 

1 + cos (vtq) 
cosh(x) + cos(7ra) 

is the characteristic function of the random variable with density 

^ (1 -|- cos(7ra)) sinh(7rat) 
sin(7ra) sinh(7rt) 

For every n > 2, si, . . . , G M and < xi < . . . < Xn, this yields 



^ - 2(1 + cos(7ra)) ^ 

l<i,j<n 



k=l 



,iykXa 



> 0, 



with the notation ti = Si/xi and yi = log(xj). Hence, the quadratic form [-/^^(xj, a;j)]i<jj<„ is 
positive definite for every n > 2 and < xi < . . . < Xn- In particular, one has 

det[Ko,{xi,Xj)]^^^ .^^ > 0, 

for every n>2 and < xi < . . . < Xn- This shows the required equivalence. 

□ 



3. Proof of part (b) 

Our argument for this part is the same as the one we have used in [13], in the framework of 
confluent hyper geometric functions. Let X = {xi < . . . < x^} and Y = {yi < . . . < be two 
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sets of positive variables and 



= n ~ 



and 



Y 



n ~ 

l<i<j<n 



be the usual Vandermonde determinants. Consider U'^{X,Y) = det^ [Ka{xi,yj)] and the afore- 
mentioned derivative determinant 



A^{x,y) = det. 



Using repeatedly the formula 



x,y > 0. 



e ^ 0+ 



fc=0 



which is valid for any smooth real function /, elementary operations on rows and columns show 
that 



(7) 



Al{x,y) = sf(n-l)^ lim 



where we have set xf = x + {i — l)e, yf = y + {i — l)e for i = 1 . . . n, = (xf , . . . , x^), Yf, = 
iVi, ■ ■ ■,yi), and 

k 

sf{k) = Ylil 
i=0 

for the superfactorial number. By Proposition 2, this entails that A^{x,x) > for any x > and 
below it will be established that the inequality is actually everywhere strict - see Remark 5. 

On the other hand, if A^(1,0+) < for some k G {2, . . . then ([7]) entails that Ka is not 
TP„ and hence not SR„ by Corollary 1. We will prove that this is the case as soon as a > 1/n and 
a ^ {1/2, 1/3, ... , 1/n}. More precisely, setting 

sin kna 



Ut 



and 



sm-TTQ 



k=l 



for every k,n > 1, we will show that 



(8) 



AS(1,0+) = sf(n-l)V-. 



We give four different proofs of each corresponding to a specific approach to evaluate the deriv- 
ative determinant A^{x,y) in closed form. The first two proofs involve classical tools, respectively 
Wronskians and Schur functions. The last two proofs rely on more elaborate and recent results 
on the Izergin-Korepin determinant, using combinatorial formula due respectively to Lascoux and 
Propp. The latter formula, involving alternating sign matrices, will be also used to obtain Part (c). 
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3.1. Wronskians. We fix a G [0, 1), set = cos(7ra) and = sin(7ra). We use the notations 

Qa{x,y) = y + xcosina), Qa{x,y) = Qaiy,x), Raix,y) = Qcix,y)Qa{x,y), 

and Pa{x, y) = x'^ + 2xy cos{Tra) + y^. In the following, we will often skip the dependence in (x, y) 
for concision. Observe that 

(9) Ra - CaPa = s^xy and Q^- P^ = -s^y'^. 

Introduce the kernels 

k=0 

and 

k=0 

for any integer p. 

Lemma 1. For every n > 0, one has 



QyU 

Proof. The formula is obviously true for n = 0, 1. By induction, it suffices to show that 

- ia,2p+l and — - Ia,2p+2 

for any integer p. The latter are elementary computations, whose details are left to the reader. 

□ 

Lemma 1 allows to express as the Wronskian of {Ta^o, ■ ■ ■ ,Ta^n-i) built on x— derivatives, in 
other words one has 



(10) AS = det. 

Observe that 



dx 



i-l 



T.,.,(1,0+) = (2p)!X:(-l)^C,V.(2cosM)^--^^ = (2p)! ""^^f + 

^ — ' Sm TTCt 

k=0 ^ ' 

where the second equality comes after identifying two standard definitions of the Chebyshev poly- 
nomial of the second kind U2p- Making the same identification for Ta,2p+i{^^^+) and putting the 
two together entail 

(11) r,„(i.o+) = (-i)V!?i*l±iM 

sm(7ra) 

for every r > 0. We next compute the successive x— derivatives of T^^r at (1,0+). The outline of 
the proof is analogous to Lemma 1, but it is more involved and so we give the details. 
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Lemma 2. For every j > 1 and r > 0, one has 

(12) = (-i)'^77TT)r^-(i'0+)- 

Proof. We first consider the case r = 2p. Setting 

_ + i + 1)! j^a _ ( — 1)-' d^Ta^2p j^a _ ( _-^\p-k (-iP-k .k r>.2k-l (-ia _ rpa 

"J'^ ~ (r + 1)! ' -^'^ ~ {2p)\ dxi ' ~ "-"p+fc^ i "-^fc.p — ^a^k,p^ 

we will show by induction that for every j > 1, one has 

(13) i?^, = aj,Mj2V'^'^'^' + yH^,P 

k=0 

for some smooth function Hj^^ on (0, oo) x [0, oo). Specifying to {x, y) = (1, 0+) clearly entails (fT2]) . 
We begin with the case j = 1, in computing 



\k=0 k=l 

V 



(2p + 2)Q„J^G^,,Kr'=+2 + IxyslY^Fl^Kl^"^^^ 

k=0 k=l 

i2p + 2)Q^j2Vc.^'^' + yHi,P 



k=0 

for some smooth function Hf^ on (0,oo) x [0, oo), where in the second equality we used the first 
identity in ([9]). Suppose now that (|13p holds for some j > 1. Differentiating, we obtain 

p 3H'^ 
= a,M-^Y.^2{p + k + j + l)Gl^Ql-2kc^P^Fl^Q^-3Gl^)Kl^^^^ - y-^ 

k=0 

= {2p + j + 2)a,M^' jZ Gt,pK^^^'^'+' + 2{R^ - cMQi E kF^,,K^^^'^^+^ 

k=0 k=0 
p f)H°' 

+ {Ql-PM'-'Y.GipK''^^'^^^" - y-^ 

k=0 



k=0 

for some smooth function -ff^+i p on (0, oo) x [0, oo), where in the second equality we used the two 
identities in ([9]). This completes the proof of (I13p for r = 2p. An analogous induction shows that 

(-1)-^' d^Ta^2p+l _ „ ^ .p+l_.fc „p-fc /.^2xfc7^p+fc+i+2 , ra 

j2pTT}\ dx^ ~ ^"j.2P+l^"^aZ^V ^) + y^j,p 

k=Q 

for every j > l,p > 0, where /"^ is a smooth function on (0, oo) x [0,oo). This yields p2p for 
r = 2p+ 1, and finishes the proof. 

□ 
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We finally deduce from (jlOp . (jlip . Lemma 2, and Leibniz's formula for determinants, that 

{i + i-l)V 



AS(1,0+) = sf(n-l)Cdet„ 



sf(n-l)2y°, 



= sf(n-l)^Cdet„^Ci,^._i 

where the last equality follows from the standard evaluation of a binomial determinant which is 
left to the reader. This completes the proof of ([8]). 



Remark 3. It does not seem that (jlOp and standard Wronskian transformations can provide any 
information which is enough explicit to characterize the everywhere positivity of (x, y) i— ?■ AJJ(x, y). 
The latter is crucial for ([2]) - see Section 4 below. 

3.2. Schur functions. This second approach relies on the following well-known expression for the 
generating function of Chebyshev polynomials of the second kind: 



\z\ < 1. 



k>0 



This entails 



n 



k>0 



with the above notation for X and Y, having written Zj = —xj^ and assuming, here and throughout, 
\yiZj\ < 1 for all i,j = 1, . . . , n. The determinant on the right-hand side, denoted by F^{X, Y), can 
be further evaluated by multilinear expansion: 

oo / n \ 

Y) = Yl det„ [4^] n K-^m' ■ 
fei...fe„=0 \/=i / 

Notice that in the multiple sum, all indices ki can be chosen distinct since otherwise the summand 

is zero. Setting 5„ for the symmetric group of size n and 

n 



1=1 



for any n— tuple K = (/ci, . . . , kn), one obtains 



FS{X,Y) = E ( E det„ [z-^'-'] (fl 

fei>->fe„>o \o-e5n \/=l 

E det„ [z';^] det, , 



(0 



ki>->k„>0 
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where in the second equahty we used twice Leibniz's formula. The above can now be rewritten in 
terms of Schur functions and Vandermonde determinants. Setting ki = n — i + Xi, one has 



n~i+Xi 
'j 



det. 



n— j+Ai 



VyVz Y.U^Sx{Y)sx{Z) 



Ai>--->A„>0 A 

where in the second equahty the sum is meant on all partitions of size n and we use the standard 
notation for Schur functions, displayed e.g. in Chapter 4 p. 124 of Putting everything together 
and differentiating with ([7]) yield after some simplifications 

sf (n - 1)2 ^ " ^ 



2.n(n+l)/2 



A \j=l / 

for all < y < X, with the notation |A| = Ai + • • • + A„. By the change of variable fii = Xn+i-i, the 
classical formula 

Ul<i<j<n(.>'^ - + j - i) 

si(n — 1) 

which can be recovered from the Jacobi-Trudi identity - see e.g. Proposition 4.2 in [3] - and a 
standard binomial determinant evaluation, entails finally 



(14) A^(x,y) 



fc=0 Ml < • ■ ■ < \«=1 / l<i<j<n 



(flj - IXi+j - if. 



fil + ■ ■ ■ + fj.,1 = k 

We can now compute explicitly AJJ(1,0+) because the summation is made on the single partition 
(0, . . . , 0) : we get 



AS(i,o+) = n X n (j'-^) 

l<j<jf<n 



sf(n-l)^y„". 



i=l 



as required by 



Remark 4. Because of the alternate signs, the above expression ()14p does not seem very helpful 
either to study the everywhere positivity of A2.{x,y). 

3.3. Rectangular matrices. This third approach hinges upon a certain closed expression for the 
determinant 

The latter is called the Izergin-Korepin determinant in the literature, and appears in the context of 
the six- vertex model - see Chapter 7 in [3] . We use its evaluation in terms of rectangular matrices 
separating the variables, which is due to Lascoux - see Theorem q in [9\. It is given by 



^^^^ det^iH^xEl^] 
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where Fq{X,Y) = ,<n(^i + 2/i)(9^i + Vj)^ and 



H 



X 



[hk-iiX)] 



l<i<n,l<fc<2n-l ' 



9^1 



l<k<2n-l,l<j<n 



are two rectangular matrices involving the complete resp. elementary symmetric functions, which 
we recall to be defined through the generating functions 



(15) 



J^hk^t" = 11(1 - Xrt)-' and J]efc(y)t'^ = 



fc>0 



r=l 



fc>0 



r=l 



Setting q = e^'™, = {xq ^^"^j . . . ,xq ^/^), and Yy = {y, . . . ,y), Lascoux's formula and ([7|) yield 

A2(x,y) = sf(n-l)2g-"("-i)/4K„(x,y)"'det„ 
By (fT5]) . we have hr{Xx) = C'!^^j._iX^q~'''^'^ , whence 



(-ik-i {i-k)/2 k-i 



l<i<n,l<k<2n-l 

On the other hand, (|15p entails CriYy) = C^y^ , so that after some simplifications 

f^n-k+j-l (k-l)/2TTa „.n~k+j~l 
y '-^k+2-2jy 



I " ' """ J l<fc<2n-l,l<j<n ' 

with our above notation for The i-th row of the product H-^q x Ey^ having a factor g'-*"^-'/^, 
we finally obtain 

(16) A^(x, y) = sf(n - 1)^ K„(x, yf det„ [A„(x) x i?-(y)] 
with the notations 

and 



An{x) 



r^k-i rrk-i 



l<j<n,l<fc<2n-l 



L " ''+^ Jl<fc<2n-l,l<i<n 

Setting L"{x,y) = det„ [A„(x) x i?"(y)] , the Cauchy-Binet formula entails 
(17) 



L^ix,y) = A.(x)B-(y) 

l<(Ti<...<(T„<2n— 1 

where A^{x) is the n x n minor obtained from the columns cJi, . . . , (T„ in An{x) and B°(y) is the 
n X n minor obtained from the rows cri, . . . , (T„ in B^{y). By Leibniz's formula, one has 

A^(x) = A^(l)x'^- and B°(y) = B°(l)y"('^-i)-"- 



with the notation 



n„ 



cr,; - I 
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This shows that L'^{x,y) is a homogeneous polynomial of degree n(n — 1) with coefficient 

(18) E A.(1)B?(1) 

for the term x^y'^^'^~^'>~^ . Besides, by (jl6p and symmetry, we have L'^{x, y) = L'^{y, x) so that these 
coefficients are palindromic. We compute 

L^{1,0+) = L^{0+,1) = Aj-,(1)B^^(1) 

with the notation Id = (n, . . . , 2n — 1). One finds immediately B^(l) = U" x • • • x and some 
elementary linear transformations yield A^(l) = 1. Since Ka{0+, 1) = 1, we finally deduce 

3.4. Alternating sign matrices. In this last approach we use Izergin-Korepin's original formula 
for the determinant Dq{X, Y) and its expression in terms of alternating sign matrices, due to Propp, 
which is given in Exercise 7.2.13 p. 244 in [3j. It reads 

(19) D,(x,y) = E(-i)-(^)(i-g)^-(^)g^-^(^)n-r^^^^f^^^ n(«^.-^+%-) 

aij = 

where An stands for the set of n x n alternating sign matrices (ASM) viz. those matrices made 
out of Os, Is and —Is for which the sum of the entries in each row and each column is 1, and the 
non-zero entries in each row and column alternate in sign. We refer to [3j for a comprehensive 
account on this topic, and also to Section 3 in the recent paper [1] for updated results. In (fT9]) . 
the following notations are used: fJ-iiA) resp. /U*(yl) is the number of —Is in the i-th row resp. i-th 
column of A, ^{A) the total number of —Is in A, I{A) the generalized inversion number of A viz. 

I {A) = ^ aijUki, 

i<k,l<j 

and 

1 otherwise. 

Notice that ASM matrices without —1 are permutation matrices. In particular, one recovers Bor- 
chardt's formula jS]) from (jl9p in setting q = 1 - see Exercise 7.2.14 in [3j. In the following we 
will set J {A) = tJ{(i,i), E k<i^kj — J2i<j'^ii}^ which is the number of southwest or northeast 
molecules in the terminology of the six-vertex model, and satisfies the formula 

(20) J{A) = 2I{A) - 2fi{A) 

(see Exercise 7.1.8 in [3J). Setting q = e^"^", we see that ([7]) entails after some simplifications 
hinging upon (|2Up the following expression for the derivative determinant A'^{x,y): 

(21) sf{n-lfK^{x,yf (P«(x, y))'^(^) (Q„(x, y))^(^) (Q„(x, y))"("-i)-2^(^) 
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with the notations 

P^{x,y) = 4sm^{TTa)xy and Qa{x,y) = e^^/^x + 6"'^/^^^ 

This yields 



A^(1,0+) = sf (n - l)2e-^""('^-i)°/2 ^ g- 



>)2p-i7rn(n-l)a/2 p2i7ra/(A) 

where 5„ stands for the set of permutation matrices of size n. Using the generating function of I{A) 
for permutation matrices given e.g. in CoroUary 3.5 of [3], further trigonometric simphfications 
entail 

A-(1,0+) = sf(n-l)V„", 

as required by ([8]). 

Remark 5. The formula (1211) also shows that 



A-(x,x) = '^^^'^ y (4sin2(7ra/2))^(^) > 

^ (2xcos(7ra/2))"('^+i) ^ v / 



for all X > and a S [0, 1). 



4. Proof of part (c) 



This last part of the theorem relies on Karlin's ETP criterion on the derivative determinant (see 
Theorem 2.2.6 in [7]) which states that if A^(x,y) > for every k G {2, . . . ,n} and x,y > 0, then 
Ka is STP„. By an induction, we hence need to show that 

(22) Q < 1/n A l/(n^ -n- 6)+ ^ AJJ(x,y) > foran2;,y>0. 

This will be obtained from (j2ip and some considerations on ASM matrices, all to be found in [3] 
and Section 2.1 of [Ij. Again, for concision we will skip the dependence of the involved kernels in 
{x,y). For any A £ An, introduce the statistics 

HA) = ^'I'^i'r = ^ 

1 < i < i'< n 
l<j<j'<n 

(see pp. 5-6 in [T] for an explanation of the second equality). Set ^„ = max{fi{A), A £ An} and 
notice that Hn = {n — 1)^/4 for n odd and that /i„ = n(n — 2)/4 for n even (see again p. 6 in [Ij). 
From (j2ip . it is clear that if 



(23) F^^, = Yl Q2KA)Q^(-i)-2KA)-2fc > 



ti{A) = k 
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for all k = . . . then is everywhere positive. Notice first that F"^. is real since it writes 

1 ^ |-Q2i/(A)Qn(n-l)-2i/(yl)-2fc _|_ Q2jy(A)Qn(ra-l)-2i/(A)-2fe-j 
^ A<£A„ 

Indeed, setting A'^ for the anticlockwise quarter-turn rotation of A, one has A'^ £ An with 
= n{A) and 2iy{A'^) = n{n - 1) - 2z^(yl) - 2^i{A) (see p. 7 in [T]). 

Suppose now that A: 7^ 0. Then necessarily n > 3, > 1 and v{A'^) > 1 (see p. 7 in fl]) and 
one gets the factorization 

ti(A) = k 
i=l 

where ^ are non-negative integer coefficients and Vn,k = max{z^(A), A G ^„,^(j4) = /c}. Notice 
in passing that no closed expression for ^ or even Un^k seem available in the literature. For n = 3, 
necessarily k = 1 = 1/3^1 and there is only one summand, so that F^^ = \Qa\^ > 0. For n > 4, it is 
sufficient to check that 

for all i = 1 . . . i^n,k- Noticing that n(n — 1) — 2k — 2(i + 1) < n? — n — Q and recalling that 

Qa{x,y) = e^W2^ + g-iW2y^ 

this inequality becomes clearly true when a < l/{v? — n — 6) after expanding the trigonometric 
polynomials, because all involved cosines are evaluated inside [0,7r/2]. 

It remains to consider the case /c = 0. Reasoning exactly as above, for every n > 2 one obtains 
F^Q > whenever a < l/n{n — 1). The following proposition yields the optimal condition a < 1/n. 
Together with the above discussion, it concludes the proof of Part (c). 



Proposition 3. For every n > 2, one has F^q > whenever a < 1/n. 

Proof. Using again the generating function of the inversion numbers of permutation matrices, we 
obtain the explicit formula 
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By an induction argument it is hence sufficient to prove tliat 

for every n > 2. Setting a = e^''^" and recalling the definition of Qa, we write 

p+q=n-l j=0 k=0 
2n-2 

Ei 2n-2-i 
CiX y 

1=0 

for some palindromic sequence {cj, i = 0, . . . ,2n — 2} which is given by 

p + q — n — 1 
j + fe = i 

i<2p 

fe<2g 

Summing appropriately and using some standard trigonometry, one can show that 

cos(7r(n — 1 — i)a) — cos(7rna) 
sm [TTa) 

for every i = 0, . . . , 2n — 2. We omit the details. This completes the proof since the latter expressions 
are all positive whenever a < 1/n. Notice also that these formulae show that the sequence {cj} is 
unimodal. 

□ 



5. Open questions 

In this section we give some heuristic reasons supporting the validity of the inclusion 
(24) a < 1/n AS(x,y) > for all x,y>0 

which, by the ETP criterion and an immediate induction, is enough to show ([2]). First, we formulate 
a conjecture on the positivity of certain generating functions of ASM matrices which would entail 
()24p . and we test it on the values n = 3,4,5. Second, we discuss more thoroughly Lascoux's 
factorization and settle two problems on the positivity on certain minors of rectangular matrices 
involving Chebyshev polynomials, whose solution would again entail ()24p . 

5.1. Alternating sign matrices. By ()2ip . the positivity of amounts to that of 
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with the notations of Section 4. The above sum can be expressed with the bivariate generating 
function 

AeAn 

which is itself given as a certain functional determinant - see Formula (57) p. 21 in [Ij. Unfortu- 
nately the latter determinant is in general very difficult to evaluate, except in certain particular 
cases. Its value at x = y = 1, counting the cardinal of An, was the matter of a whole story which 
is told in the book [3j. As in the proof of Part (c) let us now write 

PI,{A)q2u{A)q2u{AQ) ^ 
AeAn k=o 
recalling the notations Hn = max{fi{A), A G An} and 



AeA„ 

M(A) = fe 



In view of Proposition 3, it is natural to raise the following question. 

Conjecture 1. For every n >2 and every k = 0, . . . , /U„, one has F^^ > whenever a < 1/n. 
The kernel -F"^ can be expressed in terms of the generating function 



kl dy^ 

but unfortunately no tractable closed formula is known for the latter. Even its value at x = 1, 
which counts the number of elements of An with k negative entries, is known in closed form only 
in some cases - see Chapter 3 in [lOj and the references therein. 

By Proposition 3, Conjecture 1 is true for /c = 0. In the preceding section, we also proved that 
F^]^ > whenever a < 1/3. Let us now show the validity of Conjecture 1 for n > 4 and k = 
Recall that fi2 = and ^3 = 1 so that there is no loss of generality in considering n > 4. 

Proposition 4. For every n > 4 one has F"^^^ > whenever a < 1/n. 

Proof. First, suppose that n = 2p + 1 is odd. Then /i„ = and this concerns one single matrix 
with inversion number 1/ = p{p + l)/2 (see p. 6 in [Ij). We deduce 

Zn,,ni^) = and = > 0. 

Second, suppose that n = 2p is even. Then /i„ = p{p — 1) and this concerns two matrices with 
inversion numbers = p{p + l)/2 resp. = p{p — l)/2 (see again p. 6 in [IJ). We deduce 

Zn,,„{x) = x^^P~'y\l + xn and = + > if a < 1/n. 

□ 
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Since there is no closed formula for F"^. if k {0,^„}, it seems quite difficult to establish 
Conjecture 1 for all n, k. Let us conclude this paragraph in checking its validity for n = A and 
n = 5. Since fi^ = 2 and //s = 4 we just have to consider the cases k = 1 resp. k = 1,2, 3. 

• Z4,i(x) = 2x(l + xf and FJ^^ = 2|Q„|4(Q2 + q2)3 > q if a < 1/2. 

• ^5,1 (x) = x(3 + 14x + 35x^ + 48x^ + 48x*^ + 35x^ + 14x'' + 3x'^), which rewrites 

X^ — 1\ n/x^— 1\ q/x^ — 1\ A f X"^ — 1 

' " ^ / 111 n^J / I I o^4 ' 



and yields 



-^3,0 -^2,0 -^3,0 



The latter is positive for a < 1/5, by Proposition 3. 

• ^5,2 (2;) = x(2 + 12x + 21x2 + 24x^ + 21x*^ + 12x^ + 2x'^), which rewrites 

2x(l + x)(l + X + x^)(l + X + x^ + x^) + 6x^(1 + x^f + 11x^(1 + x^) 

and yields 

F5°2 = 2|Q,|Xo + Q\Qaf{Qt + Qtf + n\Q^\'\Qi + Qi). 

The latter is positive for a < 1/4, again by Proposition 3. 

• Z5,3(x) = x(l + Gx^ + 6x^ + x^) = x(l + x)(3x2 + (1 + x^)'^ - x(l - x)^), so that 

i^5"3 = \Qa\HQl + QlmQaf + {Qt + Qif-\Qa\HQl-Ql?). 

Since {Ql - Qlf = -4(x - yf sin2(7ra), we see that F^.^ > if q < 1/2. 

Observe that the generating function ^5^3 has vanishing terms so that F^^ cannot be suitably 
written in terms of F^q,F2q,F^q,F2q, contrary to F^2 -^5*1- general, testing the positivity 
of -F"^ seems to depend on bizarre rearrangements. 

5.2. Rectangular matrices. In this paragraph we study in more details the minors Ao-(l) and 
B"(l), with the notations of Section 3.3. Indeed, formulas (fT6|) . (fT7|) and ([T8|) show that the 
positivity of is ensured by that of Act(1) and B°(l). The analysis for Act(1) is easy. 

Proposition 5. For every a : {1, . . . , n} — >■ {1, . . . , 2n — 1} increasing, one has Ao-(l) = A^(l) > 0. 

Proof. The generating function 

Erir r _ ^ 

r>0 ^ ' 

and Edrei's criterion - see Theorem 1.2 p. 394 in [7j - show that the sequence {C^_x+r; r > 0} is 
TPoo. In other words, the matrix An{l) is TPqo viz. all its minors are non-negative. In particular. 
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all coefficients A„{1) are non-negative. We now compute its exact positive value with the help of 
a formula of Gessel and Viennot on binomial determinants. Substracting the i-th row from the 
(i — l)-th row successively for i = n . . .2 and then repeating this operation for i = (n — 1) ... 2, 
i = (n — 2) ... 2, ... we see indeed from the Pascal relationships Cr^^ = C^^l — Cr that the minor 
Ao-(l) is actually taken from the matrix 

^k—i 



r" 



l<i<n.l<fe<2n-l t*^^] 0<i<n-l,0<fc<2n-2 " 



The alternative formula mentioned at the bottom of p. 308 in [6j entails 

^ n.,..,,„(^-..) , 

sf (n — 1) 

In particular A„(l) = A5-(l) is always a positive integer (equal to one if and only if (7^ — o"i = n — 1). 

□ 

The analysis for B"(l) is however much more delicate. After transposition, we see that it is the 
n X n minor obtained from the columns fii, . . . , (T„ in the horizontal matrix 

Again, the generating function 

and Edrei's criterion show that the matrix [C'n~*]x<j<„ i<k<2n-i '^-^00- since some are 
negative, nothing can be said a priori about the non-negativity of B°(l). Transforming the rows 
(Li, . . . , Ln) through the simultaneaous linear operations 

n 

Li — ;> ^ C^Lj 

j=i 

multiplies the nxn minors by a constant positive factor and yields a matrix whose (i, A;) coefficient 
is given by 

E]=i Cict'U:^j_^ iil<i<n-l,i<k<n + i, 
Cn'^'U^n-k iii = n,n<k<2n-l, 

otherwise. 

Observe that the (z, n + i) coefficient equals 

n [(n+«)/2] 

E/^j /^n+i-jjja _ ST^ rij (^n+i-j /Tj-a 1 rr" ^ — n 

^n^n ^2j-n-i — ^n^n W2j-n-i "r '~'n+i-2j) — 

j=i j=i 

SO that we have a band-matrix of width n, whose (i, k) coefficient is given by 

E]=^ CiC^-^U^j_k iil<i<n,i<k<n-l + i, 
otherwise. 
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We set Ba^n for the above matrix, whose coefficients are all non-negative when a < 1/n. We also 
remark that -Ba,n is persymmetric, viz. Ba^nih k) = Ba^n{n + 1 — i,2n — h) for any (i, k). Indeed, 



a 

2n-2j-k 



— 2-^ ^ri^n '^2n-2j-k — '^n'^n ^2n-2j-k 

j=0 j=n—k 

and the last sum clearly vanishes. This persymmetry entails Bo-(l) = 6^(1). We now state a natural 
conjecture which would entail (j24 



Conjecture 2. For every n > 2, all coefficients Bct(1) are positive whenever a < 1/n. 



The difficulty to prove this conjecture comes from the rectangular shape of the matrices {B^^n}, 
which makes it seemingly impossible to use any kind of induction argument. Let us finally display 
the five ffist elements of the sequence {Ba,n}, whose common shape reminds that of Trinidad and 
Tobago's national flag (the red background being the zeroes, the crucial white diagonal stripes 
Chebyshev polynomials of the second kind, and the black central parallelogram positive linear 
combinations thereof) and displays some spatial unimodality around the middle. 



Ba,3 



Ba,l = (Ui), Ba,2 



2[/f 




u. 



2 





2C/f 




B, 



olA 



/4[/f 



V 








4^73" + 24C/f + IGU^ 




3C/2" 






U§ + Wf 3U^ 



Tja 




4C/3° + 24C/f 



3C/2" 3f/f 



4C/3" 



\ 




4c/fy 



and 



Ba' 



5U? 







V 



IOC/2" iotas' + 50t/f 5^74" + 5OC/2" 



Ui + 25Ui 









wu§ 






5Ui + 50Ui Ui + 25Ui + 100[/f 



5Ui 




u^ 



25Ui 



Ui 



Using some elementary trigonometry, for these first five values of n we could check that B, 
(viz. all its minors are non- negative) if a < 1/n The following is hence natural 

Conjecture 3. For every n > 1, the matrix Ba^n is TP whenever a < 1/n. 



\ 




5C/f/ 

a,n is TP 



This conjecture is stronger than Conjecture 2 since it involves all minors. Notice that several 
criteria have appeared in recent years to prove the total positivity of a given matrix without checking 
every minor. See all the results mentioned in Section 2.5 of [11] and also Theorem 2.16 therein for 
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a criterion only on 2 x 2 minors, interestingly also related to the zeroes of Chebyshev polynomials 
of the second kind. Unfortunately, none of this criteria seems particularly helpful in our situation. 



6. Proof of the corollary 



Setting fa{x) = Pa{^,x), we see from 1^ that pa{t,x) = t ^^"'faixt ^^'^), so that by Theorem 
1.2.1. p. 18 in jYj the kernel pa{t,x) has the same sign-regularity over (0, +oo) x (0,+oo) than the 
kernel fa{e^~y) over M x M. In Paragraph 7.12.E p. 390 of [7] it is shown that 



and the question is raised whether /o(e^~^) should be TP of some finite order when a is not the 
reciprocal of an integer. In [12], we obtained the equivalence fa{e^~^) G TP2 a < 1/2. 

Let us now prove the Corollary. We need to show that if a is not the reciprocal of an integer 
and a < 1/n, then fa{e^~^) SR„. If this were true, then the kernels /^(e^"^) and 



would also be SR„ by Theorem 1.2.1. and Lemma 3.1.1. in [7]. However, a well-known fractional 
moment identification - see (3.1) in [12j and the references therein - shows that the latter kernel 
equals Qaix — y) with the notations of Section 2. Hence, we get a contradiction to Part (b) of the 
theorem. 



We finish this paper with a natural conjecture on the total positivity of the positive stable kernel 
Pa, which reformulates Karlin's question in a more precise manner. By Lemma 3.1.1. in [7J and 
the same argument as in the proof of the corollary, this would also show the conjecture stated in 
the introduction. But we believe that this last conjecture is harder because the kernel pa is not 
explicit in general. 

Conjecture 4. For every n > 2 one has 



fa{e^~y) G STPoo ^ faie'^-y) G TPoo ^ a G {1/2, 1/3, . . . , 1/n, . . . 



} 




□ 



Pa G STP„ Pa £ SR„ 44> a G {1/2, 1/3, . . . , 1/n, . . .} or a < 1/n. 
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